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Abstract — The analysis of modulation schemes for the physical 
layer network-coded two way relaying scenario is presented 
which employs two phases: Multiple access (MA) phase and 
Broadcast (BC) phase. It was shown by Koike-Akino et. al. that 
adaptively changing the network coding map used at the relay 
according to the channel conditions greatly reduces the impact 
of multiple access interference which occurs at the relay during 
the MA phase. Depending on the signal set used at the end 
nodes, the minimum distance of the effective constellation at the 
relay becomes zero for a finite number of channel fade states 
referred as the singular fade states. The singular fade states fall 
into the following two classes: The ones which are caused due 
to channel outage and whose harmful effect cannot be mitigated 
by adaptive network coding are referred as the non-removable 
singular fade states. The ones which occur due to the choice 
of the signal set and whose harmful effects can be removed 
by a proper choice of the adaptive network coding map are 
referred as the removable singular fade states. In this paper, 
we derive an upper bound on the average end-to-end Symbol 
Error Rate (SER), with and without adaptive network coding at 
the relay, for a Rician fading scenario. It is shown that without 
adaptive network coding, at high Signal to Noise Ratio (SNR), 
the contribution to the end-to-end SER comes from the following 
error events which fall as SNR ^ : the error events associated with 
the removable singular fade states, the error events associated 
with the non-removable singular fade states and the error event 
during the BC phase. In contrast, for the adaptive network coding 
scheme, the error events associated with the removable singular 
fade states contributing to the average end-to-end SER fall as 
SNR ^ and as a result the adaptive network coding scheme 
provides a coding gain over the case when adaptive network 
coding is not used. It is shown that for a Rician fading channel, 
the error during the MA phase dominates over the error during 
the BC phase. Hence, adaptive network coding, which improves 
the performance during the MA phase provides more gain in a 
Rician fading scenario than in a Rayleigh fading scenario. Also, 
it is shown that for large Rician factors, among those removable 
singular fade states which have the same magnitude, those which 
have the least absolute value of the phase angle alone contribute 
dominantly to the end-to-end SER and it is sufficient to remove 
the effect of only such singular fade states. 

I. Background and Preliminaries 

The wireless two-way relay channel (Fig. [TJ in which 
bidirectional data transfer takes place between the nodes A and 
B with the help of the relay R is considered. All the three nodes 
are assumed to have half-duplex constraint, i.e., they cannot 
transmit and receive simultaneously in the same frequency 
band. The Denoise and Forward protocol, introduced in fl], 
is considered which consists of the following two phases: 
the multiple access (MA) phase, during which A and B 
simultaneously transmit to R and the broadcast (BC) phase 
during which R transmits to A and B. Network coding map. 



also called the denoising map, is employed at R in such a way 
that A (B) can decode the message of B (A), given that A (B) 
knows its own message. 
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(a) MA Phase (b) BC Phase 

Fig. 1. The Two Way Relay Channel 

A. Background 

The concept of physical layer network coding has attracted 
a lot of attention in recent times. The idea of physical 
layer network coding for the two way relay channel was 
first introduced in Q, where the multiple access interference 
occurring at the relay was exploited so that the communication 
between the end nodes can be done using a two phase protocol. 
Information theoretic studies for the physical layer network 
coding scenario were reported in p), H. A differential mod- 
ulation scheme with analog network coding for bi-directional 
relaying was proposed in Q. The design principles governing 
the choice of modulation schemes to be used at the nodes for 
uncoded transmission were studied in |6J. An extension for 
the case when the nodes use convolutional codes was done 
in f7|. A multi -level coding scheme for the two-way relaying 
scenario was proposed in |8|. Power allocation strategies and 
lattice based coding schemes for bi-directional relaying were 
proposed in |9|. 

Error analysis of the two-way AWGN relay channel with 
physical layer network coding based on the Detect and For- 
ward (DF) protocol, in which the relay R transmits the estimate 
of the Exclusive OR (XOR) of A's and B's transmission 
bits, was done in pO| . Performance analysis for the two-way 
Rayleigh fading relay channel with physical layer network 
coding, based on the Amplify and Forward protocol was 
presented in [ll| . For a two-way Rayleigh fading relay channel 
with BPSK modulation, upper and lower bounds on the 
Symbol Error Rate (SER) for the DF protocol were obtained 
in 1 12|. Exact BER analysis of the two-way Rayleigh fading 



relay channel with BPSK modulation for the DF protocol was 
done in |13|. 

While for BPSK modulation, the XOR network code offers 
the best performance, for other signal sets, changing the 



network coding map adaptively according to channel con- 
ditions provides a significant performance improvement |6). 
A computer search algorithm called the Closest-Neighbour 
Clustering (CNC) algorithm was proposed in |6| to obtain such 
adaptive network coding maps resulting in the best distance 
profile at R. An alternative procedure to obtain the adaptive 
network coding maps, based on the removal of deep channel 
fade conditions called singular fade states using Latin Squares 
was proposed in |14J. A quantization of the set of all possible 
channel realizations based on the network code used was 
obtained analytically in fTSl. For the adaptive network coding 
schemes, performance improvement results due to a proper 
choice of adaptive network coding maps which effectively 
mitigate the effect of most of the singular fade states, referred 
as the removable singular fade states. In |iT4J, a method 
to obtain adaptive network coding maps which remove the 
harmful effect of these removable singular fade states, using 
the completion of partially filled Latin Squares was proposed. 

UnUke the DF protocol, in which XOR map is used irre- 
spective of channel conditions, the average SER analysis of 
the adaptive network coding schemes proposed in |6J, [14j - 
ifTS) should take into account the fact that the network coding 
maps used depend on the channel fade coefficients. 

In this paper, we derive an upper bound on the end-to-end 
SER for the adaptive network coding schemes ( Q, p4| - fl5\ ) 
as well as for the case when adaptive network coding is not 
used, for a Rician fading scenario. With every singular fade 
state, we can associate an error event during the MA phase. It 
is shown that without adaptive network coding, at high Signal 
to Noise Ratio (SNR), the contribution to the end-to-end SER 
comes from the following error events which fall as SNR^ : 
the error events associated with the removable singular fade 
states, the error events associated with the non-removable 
singular fade states and the error event during the BC phase. 
In contrast, for the adaptive network coding schemes proposed 
in Q and p4) - | [T5| , the error events associated with the 
removable singular fade states fall as SNR^^ and as a result 
the adaptive network coding schemes provide a coding gain 
over the case when adaptive network coding is not used. Also, 
it is shown that for sufficiently large Rician factors, only some 
of the removable singular fade states contribute dominantly to 
the end-to-end SER and it is sufficient to remove only such 
singular fade states. 

B. Signal Model 

Let S denote the signal set of unit energy used at A and 
B, with M = 2^ points, A being a positive integer. Assume 
that A (B) wants to transmit a A-bit binary tuple to B (A). 
Let /i : F2 — >■ 5 denote the mapping from bits to complex 
symbols used at A and B. Throughout the paper all the 
fading coefficients are assumed to be Rician distributed with 
a Rician factor K. A Rician distributed random variable X 
has a scattered component and a line of sight component, i.e.. 



the transmitting nodes. Hence, without loss of generality, we 
assume 9 = 0. 

Multiple Access (MA) phase: It is assumed that the CSI 
is not available at the transmitting nodes A and B during 
the MA phase. Also, a block fading scenario is assumed. Let 
XA = /^(sa), xb = IJ-{sb) G S denote the complex symbols 
transmitted by A and B respectively, where satSb G IFj . The 
received signal at R is given by, 

Yr = Ha\/e~sXa + Hb\/e~sXb + Zr, 

where Ha and Hb are the fading coefficients associated with 
the A-R and B-R links respectively, which follow Rician 
distribution with a Rician factor K. Note that the Rician 
factor is the power ratio between the line of sight and 
scattered components. The additive noise Zr is assumed to be 
CAf{0, cr"^), where CN{Q, u^) denotes the circularly symmetric 
complex Gaussian random variable with variance a^ . The 
average energy of A's and B's transmission is equal to E^- 
Throughout, by SNR we mean the ratio %. The ratio Hb/ Ha 
denoted as z = 76-'^, where 7 g M+ and — tt < < tt, is 
referred as the fade state and for simplicity, also denoted by 

Let Sii{Ha,Hb) denote the effective constellation seen at 
the relay during the MA phase, i.e., 

Sr{Ha,Hb) = {Haxa + Hbxb\xa,xb G S} . 

Let dmin{HA,HB) denote the minimum distance between 
the points in the constellation Sr{Ha,Hb), i.e.. 



dmin{HA,HB) = min |_ff/i (x^ - x^) + Hs (x^ - a:^s) |. 



X can be written as 



pX, 



G 



where Xq is a 



circularly symmetric complex Gaussian randonm variable with 
unit variance. Since 6* is a constant, it can be cancelled out at 



(1) 



Let AS denote the difference constellation of the signal set 
S, i.e., AS = {x- x'\x, x' € S}. 

From ([T]), it is clear that there exists values of (Ha, Hb) for 
which dmin{HA,HB) = 0. Whether for a given realization 
of {Ha,Hb), dmin{HA,HB) is zero or not depends only 
on the fade state je^^ = ^f-. The values of 76-'^ for 
which dmin{HA,HB) = are of the form — ^f^, where 
Aa;^ = XA — x'j^, Axb = xb — x'g G AS and are referred to 
as the singular fade states |T4] . Note that and 00 are also 
singular fade states which occur when Hb = or Ha — 0. 
Let Ti, denote the set of all singular fade states excluding 
and 00. 

Let (.T^,^^) e 5^ denote the Maximum Likelihood (ML) 
estimate of [xatXb) at R based on the received complex 
number Yr, i.e., 

{x% Xb) = arg min \Yii - Hax'a - Hbx'^I (2) 

Broadcast (BC) phase: Depending on the value of 76^^, R 
chooses a many-to-one map Jvi'''^ : S^ — > S' , where S' is 
the signal set (of size between M and M^) used by R during 
the BC phase. The elements in S^ which are mapped on to 
the same complex number in S' by the map M.'^^^ are said to 



form a cluster. Let {Ci, C2, ■■■, Ci} denote the set of all such 
clusters. The formation of clusters for 7^^ is called clustering, 
and is denoted by C'^ . For examples of clusterings for 4-PSK 
and 8-PSK signal sets, see | [T4) . 

For a given realization of Hj^ and Hg, the choice of 
the network coding map depends only on the fade state 
76-'^ = ^, since all the distances between the points in the 
constellation Sr{Ha, Hb) can be normalized by Ha to make 
the set of all such distances a function of only 75^^. 

The received signals at A and B during the BC phase are 
respectively given by. 



Ya = H'aVe'sXr + Za, Yb = 



EgHgXl^ 



Zf_ 



(3) 



where Xji = Ai''-^{xA,XB) G S' is the complex number 
transmitted by R. The fading coefficients corresponding to the 
R-A and R-B links, denoted by i/^ and H'^ respectively are 
Rician distributed with Rician factor K. The additive noises 
Za and Zb are CAf{0, cr^). In order to ensure that A (B) is 
able to decode B's (A's) message, the clustering C''^ should 
satisfy the exclusive law |6|, i.e., 

M'''-\xa,xb)^M''-''{xa,x'b), iorxB 7^ x'b, WxAeS. 

The node A (B) can decode the message of B (A) by 
observing Ya (Yb) through ML decoding, since A (B) knows 
xa (xb) and the map J\4^'^ satisfies the exclusive law. 

Definition 1: [l4| The cluster distance between a pair of 
clusters Ci and Cj,i 7^ j, is the minimum among all the 
distances calculated between the points Haxa + Hbxb and 
Hax'a + Hbx'b e Sr{Ha,Hb), where {xa,xb) e Ci 
and {x'j^,x'g) € Cj. The minimum cluster distance of the 
clustering C'^ is the minimum among all the cluster distances. 

A clustering C^^^ is said to remove a singular fade state 
/i G H, if the minimum cluster distance dmin{C^^'^) is greater 
than zero. 

The CNC algorithm proposed in [6] obtains the map M.'''^ 
which results in the best distance profile during the MA phase 
at R, for a given 76^^. The CNC algorithm is run for all 
possible channel realizations and a partition of the set of 
all channel realizations is obtained depending on the chosen 
network coding map. For a given channel realization, the 
choice of the network coding map is indicated to A and B 
using overhead bits. 

The CNC algorithm optimizes the entire distance profile 
instead of maximizing only the minimum distance. In some 
cases, this results in the use of signal sets with a larger car- 
dinality during the BC phase and also results in an extremely 
large number of maps. For instance, for 16 QAM, the CNC 
algorithm results in more than 18,000 maps ||6l. To solve this 
problem, an algorithm called the Nearest Neighbour Clustering 
(NNC) algorithm was proposed in |6| which maximizes the 
minimum distance alone, instead of optimizing the entire 
distance profile. 

In | |T4[ , the equivalence between the network coding maps 
satisfying the exclusive law and the mathematical structure 



called Latin Squares was shown. Network coding maps which 
remove all the singular fade states were obtained by the partial 
completion of Latin Squares. 

Consider a singular fade state h = — ^f^ G % , where Ax a 
and Axb are non-zero complex numbers which belong to AS. 



Let Ax A = xa 



and Axb = xb 



Associated with 



the singular fade state h, we have the error event that the pair 
{xA, Xb) is wrongly decoded as (x^, x'b)- The CNC algorithm 
as well as the scheme proposed in |14J remove the singular 
fade state h by placing all such pairs {xa,xb) and {x'^,x'b) 
for which Ax a — xa — x'j^ and Aa;^ — xb — x'g in the same 
cluster. 

The harmful effect of the singular fade states and 00 
cannot be removed since the pairs {xa, xb) and [xa, x'g) (and 
also [xa, Xb) and (x^, xb)) result in these singular fade states 
and they cannot be placed in the same cluster without violating 
the exclusive law |6|. The singular fade states and 00 which 
occur due to channel outage, irrespective of the signal set used, 
are referred as the non-removable singular fade states. The rest 
of the singular fade states which depend on the signal set used 
are referred as the removable singular fade states. Henceforth, 
unless explicitly mentioned, by singular fade state, we refer to 
only the removable ones. 

Throughout the paper, in a statement if it is mentioned 
simply as adaptive network coding, it refers to the scheme 
proposed in |61 as well the one proposed in fTA\, i.e., the 
claim made in the statement holds for both the schemes. 

In this paper, an upper bound on the average end-to-end 
SER is obtained for the two-way relaying scenario with and 
without adaptive network coding. From the obtained analysis, 
the reason why adaptive network coding provides performance 
improvement becomes very clear Also, it is shown that not 
all the singular fade states contribute equally to the end- 
to-end error probability. Removing only those singular fade 
states which contribute dominantly to the end-to-end error 
probability reduces the system complexity at the relay as well 
as the number of overhead bits required for indicating the 
choice of the adaptive network coding map. 

The contributions and organization of this paper are as 
follows: 

• An upper bound on the average end-to-end SER for 
the wireless two-way relaying scenario with and without 
adaptive network coding is derived (Section II A and 
Section II B). It is shown that without adaptive network 
coding, at high SNR, the contribution to the average 
end-to-end SER comes from the following terms which 
decrease as SNR^^: the error events associated with the 
removable and the non-removable singular fade states 
and the error event during the BC phase. In contrast, 
for the adaptive network coding schemes proposed in |6| 
and [14 1, the error events associated with the removable 
singular fade states fall as SNR^^ and as a result the 
adaptive network coding schemes provide a coding gain 
over the case when adaptive network coding is not used. 

• It is shown that in a Rician fading scenario, the error 
during the MA phase dominates over the error during the 



BC phase. Hence, the adaptive network coding schemes, 
which improve the performance during the MA phase, 
provides more gain in a Rician fading scenario than in a 
Rayleigh fading scenario (Section II). 

• It is shown that in a Rician fading scenario, the removal 
of the singular fade state s = 1 assumes greatest sig- 
nificance. While in a Rayleigh fading scenario, all the 
singular fade states contribute dominantly to the overall 
average SER, it is shown that in a Rician fading scenario, 
for sufficiently large Rician factors, among those singular 
fade states which have the same magnitude, only those 
for which the absolute value of the phase angle is the 
least contribute dominantly to the end-to-end SER and 
it is sufficient to remove only such singular fade states 
(Section IE). 

• Simulation results which confirm the above mentioned 
facts are presented in Section IV. 

Notations: Q[.] denotes the tail probability of the standard 
Normal distribution. Ph{E} denotes the probability of the 
event E conditioning on the set of random variables H. ]E(X) 
denotes the expectation of X. 

II. ERROR ANALYSIS OF THE WIRELESS 
TWO-WAY RELAYING SCENARIO 

In the section, upper bounds on the end-to-end SER are 
obtained for the wireless two-way relaying scenario with and 
without adaptive network coding. 

Let x^ and x^denote the messages decoded by A and 
B respectively at the end of the BC phase. Let Ea and 
Eb respectively denote the error events x^ 7^ xb and 
x^ ^ XA- Let H = {HA,HB,Hj^,Hg) denote a particular 
realization of the channel fade coefficients. The end-to-end 
SER Ph{Ea U Eb} given in Q, can be upper-bounded as in 
(|5]l (both shown at the top of the next page). 

The first and second terms in dSll, denoted as Prr' and 
Pjj' , are the probability of error events at node A and B 
respectively at the end of the BC phase, given that the relay 
decoded to the correct cluster during the MA phase The third 
term in (|5]l, denoted as P^^^, is the probability that the relay 
decodes to the wrong cluster during the MA Phase. 

The choice of the signal set S' used at R during BC 
phase depends on the number of clusters in the clustering 
C'^ f61, P41. Let Cmax denote the clustering which has the 
maximum number of clusters over all possible je^^ and let 
S' max denote the signal set associated with the clustering 
Cmax. Let P^'B^ = E(P^'^^) and P^'^^ = i?(P^'^^). 
The upper-bound on the average error probability during the 
BC phase calculated for the case when R always uses the 
signal set S' max irrespective of 76-'^ will serve as an upper- 
bound on p^.^c and p^^^c as well. Hence, we have. 



:,A,BC 



pB,BC ^^-K_ 



1 



1 

SNRd^.^(5'„ax) ' 



where dmin{S' max) is the minimum distance of the signal set 

'-^ max ■ 



The upper-bound on p^.^C and p^^^'^ given above is 
not tight. Since the performance advantage due to adaptive 
network coding is captured only by the term P^^ in IS), 
in the rest of the paper we do not focus on the probabilities 

pA,BC ^nd PB,BC^ 

The probability that the relay decodes to the wrong cluster 
during the MA Phase P^^^ can be upper-bounded as. 






^ Ph{{x%x"b)^{x'a,x'b)}. 



(x'j^,x'g):M-<-'>(x'j^,x'g)^M'<-0{xA,XB) 



(6) 



The probability Ph{{x^tX^) — {x'a,x'b)} can be upper- 
bounded by the corresponding pair-wise error probability 
given by. 



(a^A.^s)} = 

'y/SNR\HA{xA 



Q 



■ x'j^) + Hb{xb - x'l 



V2 



Hence, from (|6]l, we get ^ (shown at the top of the next 
page), where l{c} is the indicator function which is one if the 
condition c is satisfied and is zero if it is not satisfied. 

Averaging with respect to the fade coefficients in (|7|, 
we get ([8]) (shown at the top of the next page), where 
TZ{{xa,xb), {x'a,x'b)} C C^ is the region 

{{Ha,Hb):M"-^''Hxa,xb)^M"-/"Hx'a,x'b)}- 

Let P^^^{{xa,xb) — > {x'j^,x'g)} denote the term inside 
the summation in (|8]l. In the following subsections upper 
bounds on P^'^"^{{xa,xb) — > {x'a,x'b)} are obtained for the 
two-way relaying scenarios without and with adaptive network 
coding. 

A. Two-way Relaying without Adaptive Network Coding 

Consider the situation where R uses the same clustering C 
which does not remove any of the singular fade states, for all 
values of 76-'^. Since R uses the same clustering for all 76-'^, 
the region TZ{{xa, xb), {x'a, x'b)} can be either null set or the 
entire C^ plane. If {xa,xb) and {x'j^,x'b) are placed in the 
same cluster by the clustering C, then the region is the null set 
and P^^^{{xa,xb) -^ {x'a^x'b)} = 0. But since C does not 
remove any of the singular fade states, for every singular fade 
state — ^f^, Axyi, Axb € AiS, there exists at least one pair of 
two-tuples {{xa,xb), {x'a, x'b)} which satisfies Ax a — xa — 
x'a and Ax b = xb—x'b and for which {xa,xb) and {x'a, x'b) 
are not placed in the same cluster. For such pairs {xa,Xb) 
and {x'j^,x'g), TZ{{xa,xb), {x'a, x'b)} is the entire C^ plane. 
For this case, the probability P^^^{{xa,xb) -^ {x'a,^'b)} 
defined in (|6| is given in (|9]l (shown at the top of the next 
page). In (|9]l, the suffix FNC indicates that fixed network 
coding is used at R, irrespective of the channel conditions. 
Also, in Q, /{Ha) and /{Hb) denote the probability density 
functions of the random variables Ha and Hb respectively. 
Substituting Rician probability density functions for /{Ha) 



PH{EAUEB} = PH{EAUEB\M''-''{x'l,x'^)=M''-''ixA,XB)}PH{M-''\x'^,X^)=M-''''{xA,XB)} 

+ Ph{EaU EB\M-'-''ix2,S:B) ^ M''-''ixA,XB)}PH{M-'-''{x^,x'^) ^ M-'-''{xa,xb)} (4) 
< Ph{Ea\M-''''{x'^,x'^) = M'''''{xa.xb)} + Ph{Eb\M^'''{xa,Sb) = M^'''{xa,xb)} + Ph{M-'-''{x'^,x'^)^M-'-''{xa.xb)}. 



(5) 



E 



(:^',^')-.M-'.0(^',:^')^M-I'0(^^,^b) 



VSNR\Ha{xa - x'a) + Hb(xb - x'b)\ 



V2 



E 



(x' ,1' ):(x' ,x' )^(x_4,x_g) 



■{M''.0{:o' :c' J^M-I.OI^ca.^b)} 



VSNR\Ha(xa - x'a) + Hb(xb - x'b)\ 



V2 



(V) 



-jMA A -yn,/ T-.M A 



hph^) < E 



{■X ^,x ^):{x j^,x ^)j^{x j:^,x q) 



n., 



(aiA.^flJ.C^'A'^B)} 



VSNR.\Ha{xa -x'a) + Hb(xb -x'b)\ 
72 



f(HA)f{HB)dHAdHB . (8) 



pAf-4{(xA,XB)^(x' ,x' )} 



PFNci^^A,XB) -^ (x'a^^'b)} =2'^ 



^SNR\Ha(xa - x'a) + Hb(xb - x'b)\ 
72 



f{HA)fiHB)dHAdHB 



(9) 



and /{Hb), at high SNR, the integral given in (|9| can be 
upper bounded as, 



Consider error events {xa,xb) 
which XA — x'a or xb = x'b- 



— > (x'j^jx'g) at R for 
For such error events. 



P^l^c{i^A,XB)^{x'A,x'g)}< 



lAx^l^ + IAxs 



&NR{\Axa\^ + \Axb\'^) 



1 + 



(10) 



From ( (T0| , it can be seen that P^^^{(a;yi,a;B) ^■ 
(a:'^,x'^)} decreases as SNR^^ at high SNR. 
The probabilities p-4,BC ^jj^j pB.BC ^g proportional to 



the region 7?.{(a;^,a;B), (x^,a;'g)} is the entire C^ since the 
pairs {xatXb) and [x'j^^xb) (and also the pairs {xatXb) 
and [xatx'^)) cannot be placed in the same cluster without 
violating the exclusive law. At high SNR, the probabili- 
ties P^'^^{{xa,xb) -^ ixA,x'B)} and P^'^{{xa,xb) -> 
(x^, Xb)} defined in ^ can be upper-bounded as. 



-K 



-K 



at high SNR. If Axs belongs to A5, -Axb also PANc{i^A,XB) ^ [xa.x'b)} < 



belongs to A5. For Aa;^ 7^ and IS.xb 7^ 0, one of the 
two |Aa;^ + Aa;^^ and |Axa — AxbP, has to be less than 
lAx^P + lAxsp. Hence, from ^, P^'j^^{{xa.xb) ^ 



pMA 

{x'j^,x'^)} is proportional to e^'^^ where k < 1, for some 
xa 7^ x'j^ and xb 7^ x'g . When the Rician factor K increases, 
the contribution of the error during the BC phase to the 
overall average SER decreases and for large values of K, the 
contribution to the overall average SER comes totally from 
the error during the MA phase. The reason for this is that 
Pbc snd Pbc (proportional to e~^) decrease faster with 
K than those terms P'^^^{{xa,xb) — > {x'a, x'b)} which are 
proportional to e^'^^ where k < I. For this reason, adaptive 
network coding, which improves the performance during the 
MA phase, provides more gain in Rician fading scenario than 
in a Rayleigh fading scenario, consistent with the simulation 
results in [6| and |15]. The exact reason why adaptive network 
coding improves the performance during the MA phase is 
described in the next subsection. 

B. Two-way Relaying with Adaptive Network Coding 

Consider the adaptive network coding schemes proposed in 
||6) and | [T4) . In both the schemes all the singular fade states 
are removed by a proper choice of the clustering. 
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Panc{(^a,xb) -^ {x'a,xb)} < 



1 + 



SNRIAzaP 



(11) 



(12) 



The suffix ANC in the above two equations indicates that 
adaptive network coding is used at R. 

Let s — ~^^,Axa,^xb 7^ e AS denote a singular 
fade state. For the adaptive network coding schemes proposed 
in |J6) and | |l5| , the complex fade state i'^e^^) plane can be 
quantized into different regions depending on the clustering 
used at R. In the neighbourhood of every singular fade state, 
an associated region Fg exists in which a clustering which 
removes that singular fade state is used at R. For example, 
for the case when 4-PSK signal set is used at A and B, the 
quantization of the complex fade state plane is as shown in 
Fig. [2] along with the regions Ri — R12 associated with the 
12 singular fade states LISJ. The regions R13 and R14 in Fig. 
I2] are the clustering independent regions in which choice of 
the clustering does not matter and any clustering satisfying 
the exclusive law gives the same performance. For details, see 

Note 1: The quantization of the complex fade state plane 
is the same for 4-PSK signal set for the Nearest Neighbour 



Clustering algorithm proposed in |14| and for the scheme 
proposed in |6|, while it need not be the same for other signal 
sets. Nevertheless, for both the schemes, there exists a region 
in the neighbourhood of every singular fade state in which a 
clustering which removes that singular fade state will be used. 
Let Cs denote the circle with the largest radius Sg enclosed 
in the region Fj, with center at the singular fade state s. For 
example, for 4-PSK signal set, s = 1 + j is a singular fade 
state and the circle Ci+j of radius (5i+j enclosed in the region 
Fi+j (the region R3) is as shown in Fig. [2] 




Fig. 2. The Diagram showing the fade state quantization for 4-PSK signal 
set 



Let £{Cs) denote the region exterior to the cir- 
cle Cs in the complex plane. If 76^^ lies inside the 
circle Cg, M'^'''{xa,xb) = M'^-'^{x'j^,x'g). Hence, 
n{{xA,XB),{x'A,x'B)} C {{Ha.Hb) : ff e £{Cs)} = 
{{Ha,Hb) : Iff + fit I > Ss}. Hence, for the adaptive 
network coding schemes, the probability P^^^{{xa,xb) ^• 
{x'j^,x'g)} defined in ^ can be upper-bounded as given in 
( [T3| l (the suffix ANC in ( [T3] l indicates that adaptive network 
coding is used at R). Let Axa = xa — x'j^ and Aa;^ — 
xb - x'g. At high SNR, Pfj^c{i^A,XB) ~> {x'^^x'^)} can 
be upper-bounded as stated in the following theorem. Note 
that the bound is valid only at high SNR and is obtained by 
upper-bounding the term on the right hand side of ( [T3| ). 

Theorem 1: For xa 7^ x'j^ and xb 7^ x'g, at high SNR 
PAN^c{i^A^ Xb) -^ (a;^, x'^)} can be upper-bounded as given 
in ( fT4l ) (shown at the top of the next page). 

Proof: See Appendix. ■ 

From Theorem 1, it can be seen that for xa 7^ x'^ and xb 7^ 
x'g, at high SNR, P^j^(j{{x Ai x b) — >■ (2:^,3;^)} decreases 
as SNR^^. But the overall diversity order of the end-to-end 
SER will be one, since P^^Kx^i, a;^) -^ [xa^x'^)} and 
PancH^a^xb) -^ {x'j^^x'g)} given in ([TT]) and ([T2| as well 



as the probabilities p^.^C and pB,BC decrease as SNR^^ at 
high SNR. 

Even though adaptive network coding does not provide 
any diversity advantage, it provides coding gain advantage 
over the case when adaptive network coding is not used. 
For the case when adaptive network coding is not used, 
from ( [Tol l, it can be seen that Pp^j^(j{{xa,xb) -^ (a;^,x^)} 
decreases as SNR~^ at high SNR, for all pairs {xa.xb) and 
(a;^, x'g). In contrast for the adaptive network coding scheme, 
only those probabilities P^^j^fj{{xATXB) — > {^'a^^'b)} for 
which XA — x'j^ or xb ~ x'g decrease as SNR~^ and the 
rest decrease as SNR^^. In other words, for the adaptive 
network coding scheme, at high SNR, the pair-wise error 
events associated with the non-removable singular fade states 
and 00 decrease as SNR^^ and the pair- wise error events 
associated with the removable singular fade states decrease as 
SNR^^. The adaptive network coding scheme provides coding 
gain advantage over the case when adaptive network coding 
is not used, by making the probabiUty of the error events 
{xa,xb) — > {x'^,x'g) for which xa 7^ x'^ and xb 7^ x'g 
decrease as SNR"^ instead of SNR"\ 

From ( [TOJ i, for the case when adaptive network coding is 
not used, it follows that Pp^j^fj{{xA,XB) — >■ (a^^,a;^)} is 



proportional to e i+i^i^ . This suggests that removing certain 
singular fade states assumes more significance than removing 
the others and is discussed in the following section. 

III. Partial Removal of Singular Fade States 

In this section, it is shown that only some of the removable 
singular fade states contribute dominantly to the overall SER 
and only such singular fade states can be removed without a 
significant degradation in performance. Note that whether a 
singular fade state contributes significantly to the error prob- 
ability depends on the Rician factor as well. For a Rayleigh 
fading scenario (Rician factor K=0), contributions of all the 
singular fade states to the overall SER are significant. 

{{xa,xb) — > (a;^,a;^)} is proportional to 

, the factor /(s) - i^^^ 



Since Pppf(j-\ 



i+kr 



-, referred to as the domi- 



e 1+ 

nance factor of the singular fade state s, determines whether 
the contribution from a singular fade state s is significant or 
not. The lesser the value of the dominance factor, the more 
the contribution of the singular fade state s to the SER. Since 
/(s) = if and only if s = 1, the removal of the singular 
fade state s — 1 assumes greatest significance. Note that for 
the case when A and B use the same signal set, s = 1 will 
always be a singular fade state. Let ^s denote the phase angle 
ofg. We have, /(.) = !- ^l-l-|f) . 

Among those singular fade states which have the same 
absolute value |s|, those which have a lesser value of |Zs| 
(|Zs| ranges from to tt) have a lesser dominance factor 
and hence contribute more towards the overall SER. For 
sufficiently large Rician factors, among those singular fade 
states which lie on the same circle, R can choose to remove 
only those singular fade states which have the least absolute 
value of the phase angle. 



J J l(H ,.H„):\^ + ^i^^\>6.\. 



VSNR\Ha{xa - x'^) + Hb{xb - x'b)\ 



■~/2 



f{HA)f{HB)dHAdHB. (13) 



Panc{^^a,xb) ->■ (x'a,x'b)} < 



cxp •; -2K - 



A-(A- + l)|l 
|Ax_b| 



^i-<--)(-^) 



[K + 1 + SNRi- 



•)([K+l + SNR^^]5? + (K + l)(l+^f^)) 



(14) 



Example 1: For the case when 4-PSK signal set is used 
at the nodes A and B, the twelve singular fade states are as 
shown in Fig. [3] For 4-PSK signal set, the dominance factor 
f{s) for the 12 singular fade states are given by, 



s 


/(s) 


1 





J, -1,-j 


1 


0.5 + 0.5j, 0.5 - 0.5j, 1 + J, 1 - J 


1/3 


-0.5 + 0.5j, -0.5 - 0.5j, -l+j,-l- J 


5/3 



Among all the singular fade states, the singular fade state 
s = 1 is the most dominant one. Among those singular fade 
states which lie on the circle with radius 4=, the singular 
fade states 0.5 + 0.5 j and 0.5 — 0.5 j' are the dominant ones. 
Similarly, among those singular fade states which lie on the 
circle with radius \/2, the singular fade states 1 + j and 1 — j 
are the dominant ones. The relay can choose to remove only 
the dominant singular fade states on each circle, which are the 
circled ones shown in Fig. [3] 



IV. Simulation Results 

The simulation results presented are for the case when the 
end nodes use 4-PSK signal set. The fading coefficient Ha, 
Hb, H'j^, H'q are Rician distributed with a Rician factor 
K — A and unit variance. For comparison, we consider the 
case when R uses the Modulo-4 Latin Square shown in Fig. 
[5] irrespective of the channel condition (every entry of the 
Modulo-4 Latin Square is the modulo 4 addition of the row 
index and the column index). Note that the Modulo-4 Latin 
Square does not remove any of the 12 singular fade states. Fig. 
|4] shows the SNR vs BER plots for the different cases. From 
Fig. p] it can be seen that the diversity order for all the cases 
considered is one. Also, it can be seen that at high SNR, the 
adaptive network coding scheme based on the removal of all 




Fig. 3. The Diagram showing the dominant singular fade states for 4-PSK 
signal set 



the singular fade states using Latin Squares proposed in 1 14| 
provides nearly 8 dB gain over the case when Modulo-4 Latin 
Square is used irrespective of channel conditions. For details 
regarding the Latin Squares which remove the singular fade 
states, see |14]. 

Note 2: For 4-PSK signal set, the adaptive network coding 
scheme based on the removal of singular fade states proposed 
in p4) and the one based on the Nearest Neighbour Clustering 
algorithm proposed in |6| turn out to be the same. 

Fig. |4] also shows the SNR vs BER plots for the case when 
only the singular fade state s = 1 is removed and for the 
case when all the singular fade states other than s = 1 are 
removed. Note that the bit-wise XOR network code removes 
the singular fade state s = 1 p4) . It can be seen that the SNR 
vs BER performance for both these cases are nearly the same. 
This means that the performance improvement provided by 
removing all the 11 singular fades states other than s — 1, 
is equal to the performance improvement provided by the 
removal of s = 1 alone. This confirms the assertion made 
earlier in Section III that the removal of the singular fade 
state s — 1 is of greatest significance. Fig. |4] shows the SNR 
vs BER plot for the case when only the five dominant singular 
fade states (the circled ones shown in Fig. l3]l are removed. It 
can be seen from Fig. HI that at high SNR, removing only 
the five dominant singular fade states results in nearly 7 dB 
performance improvement over the case when Modulo-4 Latin 
Square is used irrespective of channel conditions. In other 
words, at high SNR, the degradation in performance that 
results because of removing not all but only the dominant 
singular fade states is less than 1 dB. 
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Fig. 4. SNR vs BER plots for the different schemes for Rician factor K = A 
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Fig. 5. The Modulo-4 Latin Square 



V. Discussion 



An upper bound on the average end-to-end symbol error 
probability was obtained for the two-way relaying scenarios 
with and without adaptive network coding. From the analysis, 
the reason why adaptive network coding schemes provide 
performance improvement becomes clear Also, it was shown 
that in a Rician fading scenario, some of the singular fade 
states contribute more to the average symbol error probability. 
Simulation results show that removing only such dominant 
singular fade states results in almost the same performance as 
that of the case when all the singular fade states were removed. 
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Appendix 

Before we prove Theorem 1, we prove the following lemma. 

Lemma 1: For cq, r G M+ and he G C, the integral 

J {/ieC:|/i|>Co} '' 

where Qi is the first order Marcum Q function. Also / can 
be upper bounded as, 

1 



/ < 



Co 



-r(co-\K\Y 



rco-|/ic| 

Proof: Let 9 denote the phase angle of h and 7 denote 
the absolute value of h. 

Then the term inside the integral / can be written as, 

g-r(7" + |/ie|")g27r|h<,| cos(e-0) 

where (j) — tan^^ ( Reih ^ ) • Transforming h into polar 
coordinates (7,6'), the integral / can be written as, 

1 



/ = / 76 

7=co 



-'■(7^+l''c|"). 
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Let /n denote the Bessel function of zeroth order of the first (|24l) can be written as 
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The above inequality follows from the fact that Qi{a,l3) < 
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Using the transformation x = v2r7, the integral / given in 



PROOF OF THEOREM 1 

Consider the integral I{Ha) given in ([TSj (the equations 
([T5|l-(|23|l are shown at the top of the previous page). Substi- 
tuting the Rician probabiHty density function fiHs) in (fTSl) 

^2 — 

and upper-bounding Q{x) by e~^ , we get ( [16| ). Using the 
transformation H'g ^ Hb + Ha ^f^ in ^^^ we get (il) , 

where c^Ha^''^ ' ^ 



AxB ' V K+1- 

Using Lemma 1, the integral in ([TtJi can be upper-bounded 
as given in ([TSj. At high SNR, ([T8|, can be approximated as 
in ([191). 

The upper bound on the probability 

PANci-'-A, xb) -^ {x'j^,x'g) given in (\3\ can be written as, 

P^'j^ci^A, Xb) ^ [x'a, x'b) = / I{HA)I{HA)dHA. 

JHAec 

(25) 

Substituting the Rician pdf /{Ha), using ([T9|, ( |25] l can be 
upper bounded as given in (pOj. Rearranging the terms in ( pO| 
one gets ( pT) , where s = "sf^- Using Lemma 1, the integral 
in ( |2T] i can be evaluated as given in ( |22] i. Since (5i(a, 0) = 1 
(|22]i can be simplified as in (|23|l. This completes the proof. 



